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LIST OF SYMBOLS 



b 


thickness of a sublayer of a "slab" whose overall 
thickness is h (h < b) 


h 


thickness of the translating slab (containing par- 
ticles) through which a molecule passes (corres- 
ponds to ring thickness in application of the 
collision theory) 


n (r) 


number density of particles in a slab (or in the 
ring for applications of the collision theory) 


r 


particle radius (sometimes referred to in this 
paper as the "size" of a particle) 


-+■ 

V 

ej 


ejection velocity of a molecule, relative to its 
source particle 


V 

m 

-»■ 

V 

P 


velocity of a molecule 

velocity with which the slab translates (equivalent 
to a representative) velocity of the particles 
"contained" in the slab 


D i 


distance that the slab translates in time t 


L 


scale length for the extent of the slab in direc- 
tions normal to its thickness 


N 


equilibrium number of molecules (of a given type) 
in ring atmosphere deriving from a specific produc- 
tion mechanism 


N (r i ) 


number of slab (ring) particles of size r-j_, con- 
tained in the slab (meaningful only for discretely 
distributed r) 


P 

c 


probability that a molecule collides with at least 
one ring particle during a given passage through 
the slab (ring) 


P (r . ) 
no 1 


probability that a molecule passes through the 
slab (ring) without colliding with a particle of 
size r^ (only strictly meaningful for discrete r) 


no (3 ) 


probability that j particles of size r. are not in 
the collision region 
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Q e (r,A) 


extinction cross section, sometimes referred to as 
the extinction coefficient (relates to the effec- 
tive size of a particle with regard to removal of 
electromagnetic energy from a beam via scattering 
and absorption) 


R 


orbital radius (see Section III) 


R I 


interaction radius (see Appendix B) 


S 


sticking coefficent (the probability that a mole- 
cule will adhere to a surface when it collides 
with that surface) 


T 


a "mean" orbital period typical of the trajectories 
of all molecules of a certain type produced by a 
given process 


T 

P 


number of completed orbital periods 


W 


production rate (of molecules) associated with a 
specific mechanism, usually including only those 
molecules which assume closed orbits 


a, 8 ,y 


direction cosines which specify a molecule's path 
as it passes through the slab, when working in 
a frame of reference in which the slab translates 

at v (see Figure 2) 

P 


a' , B ' , Y ' 


direction cosines which specify a molecule's path 
through the slab, when working in a frame of refer- 
ence in which the slab is fixed (see Figure C2 ) 


4> 


angle used to^describe the deviation of the actual 
direction of v from the assumption that v is 
constant ^ ^ 


A 


wavelength (of electromagnetic radiation) 


y 


angle that an electromagnetic beam makes with the 
surface normal of a slab (ring) 


9 


angle in the x-y plane between the projection of 

the molecular path and the +x direction (Vp 

direction) , defined in a frame of reference in 
-> 

which v p 0, referred to in this paper as the 

"orbital aximuth" 


9 ' 


same as 9., except defined in the frame of reference 

with v = 0 
P 
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INTRODUCTION 



I . 

It has been well established that the particles of which 
Saturn's rings are composed are primarily water ice (Cuzzi 
et al . , 1984). Accepting this, one would expect that various 
erosive processes would produce a ring atmosphere comprised 
of H 2 O molecules and their constituents, e.g,, H atoms. The 
density of such an atmosphere would be determined by balancing 
the production rates (of atoms and molecules) associated with 
the erosion mechanisms with losses resulting from capture by 
Saturn, escape from the Saturnian system, and recapture by the 
ring particles. 

Early theoretical estimates of a ring atmosphere (Denne- 
feld, 1974; Blamont , 1974) considered sources such as subli- 
mation, meteroid impact, and bombardment by solar and 
interstellar wind.'*' These analyses were concerned only with 
the atmosphere associated with the A and B rings (see Figure 1) , 
by far the most significant potential sources of material for 
a ring atmosphere, owing to their large masses (Appendix A 
gives a summary of ring nomenclature and dimensions). Denne- 

feld calculated the density of the f^O atmosphere to be ~25 
-3 

cm , the molecules confined to a toroidal region closely 
surrounding the A and B rings. He estimated the H atom 



Dennefeld presented a conceptual model for the study of 
ring atmospheres, which is of fundamental importance. This 
model is discussed later in this section. 
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to 



Figure 1. Nomenclature of the Saturnian Ring System 
(drawn to scale) 

component of the atmosphere to have a density of ~1 cm ^ 
distributed through a sphere of ~3 xlO"* km, well beyond the 
outside edge of the A ring. 

Schemes for detection of the atmosphere all rely on 
measurement of H Lyman-a emission from the vicinity of the 
rings, the most sensitive means available. Thus, analysis is 
limited to the neutral H atmosphere. Using this approach, 
several researchers put Dennefeld's estimates to the test. 
Observations from instruments onboard rockets (Weiser et al . , 
1977), satellites (Barker, et al . , 1980; Clarke, 1981), 
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Pioneer II (Judge et al., 1980), and Voyager I (Broadfoot et 
al., 1981) detected H atoms near the A and B rings. Results 
were reasonably consistent, suggesting a H density 2-3 orders 
of magnitude greater than Dennefeld predicted. The shape 
of the atmosphere (spherical versus toroidal) cannot be con- 
clusively deduced from these Lyman-a measurements, but pre- 
liminary analysis of Voyager data indicates that it does not 
extend beyond the outer edge of the A ring (Cuzzi et al . , 

1984) . 

Efforts to reconcile theoretical estimates of the H 

atmosphere with observations have primarily been directed 

toward proposing new sources or re-evaluating those previously 

considered. Proton flux from Saturn's ionosphere (Ip, 1978), 

photodissociation of t^O on the ring particles (Carlson, 1980), 

magnetospheric ion sputtering (Cheng et al . , 1978,1980), and 

meteroid impact (Morfill et al . , 1983; Ip, 1983) are among 

2 

the most seriously studied erosive mechanisms. Of these 
meteroid impact appears likely to be the dominant source. 

While understanding of the ring atmosphere has increased, 
Dennefeld' s conceptual model has remained the foundation for 
its theoretical analysis. It is worthwhile to review this 
model. Suppose that a particular erosive process produces 
atmospheric molecules (e.g., H^O or H) with these molecules 



Some have also considered H atom sources external to the 
rings such as atoms escaping from Saturn's atmosphere 
(Shemansky and Smith, 1982). Preliminary calculations indi- 
cate that this process may be quite significant. 



11 



having some initial velocity distribution. The trajectory 
of each molecule is primarily governed by its gravitational 
interaction with Saturn, the dominant central body in the 
Saturnian system. If the ring atmosphere is quite tenuous, 
collisions between molecules as they move in their Keplerian 
orbits is unlikely (this assumption poses no problem for 

_3 

Saturn's rings, since for a density of even 1000 cm the 

mean free path is ~10^ km). A molecule's initial velocity 

and position within the ring determines its orbit for which 

3 

there are three possible cases : 

(1) the trajectory takes the molecule into Saturn's 
atmosphere, thus removing it from the ring atmosphere. 

(2) . the orbit is open (hyperbolic) with the molecule 

escaping the entire Saturnian system. 

(3) the orbit is elliptic with the molecule passing 
through the ring once or twice during each orbital 
period (depending on the eccentricity of the orbit) . 

For those molecules with elliptic orbits, some will be lost 

from the atmosphere to the rings. The probability that this 

occurs during each passage of a molecule through the rings 

is the product of the probabilities of its colliding with 

ring particles and of its sticking to a particle during this 

process. The spatial extent of the atmosphere tends to be 

toroidal if the initial speed of the molecules relative to 

the orbital speed of the ring particles is small. On the 



3 

The first two cases apply if the erosive mechanism 
produces molecules with high speeds relative to the ring 
particles; the third is dominant for relatively lower energy 
processes . 
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other hand if the initial speed is much greater than that of 
the ring particles, then the atmosphere will tend to be more 
spherical. A mathematical expression of this process of the 
form introduced by Dennefeld is 



N 




t 



where : 



N = number of molecules (of a given "kind") in 
the equilibrium atmosphere (from a specific 
production mechanism) . 

W = production rate (sec ^) of molecules asso- 
ciated with a specific mechanism which 
assume elliptic orbits (call this the 
"source term"). 

T = mean orbital period (sec) ; this term incor- 
porates initial velocity information, so it 
might more generally be referred to as the 
"ballistic term. " 

S = sticking coefficient, the probability that 
a molecule will adhere to a ring particle 
when it collides with it. 

P = probability that the molecule collides with 
a ring particle as it passes through the 

ring. ^ 

This model is admittedly crude, but it does contain the 
essence of the physical system. To see how the theory has 
evolved, each term is discussed below. 



4 

A tacit assumption in this formulation is that the 
average trajectory intersects the ring only once during 
an orbital period. If it passes through the ring twice 
during a period, P c hasa slightly different interpretation, 
but the other terms are unaffected. 
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The source term, W, has received much attention as dis- 
cussed in some detail earlier. 

The ballistic term, T, is receiving increased emphasis. 
Ballistic transport models (Durisen, 1984; Ip, 1983) incor- 
porate inhomogeneities found in the macrostructure of the 
rings with improved estimates of initial velocity distribu- 
tions to allow better specification of molecular trajectories. 

The sticking coefficient, S, is probably the least under- 
stood factor in the theory. The physical chemistry under- 
lying the phenomenon is not well developed, while laboratory 
experiments • do not appear to accurately replicate conditions 
found, in the ring environment. ^ While it may be possible that 
the physical chemistry problem will be solved, a more basic 
dilemma will likely remain for quite some time: S cannot be 

determined accurately unless detailed knowledge of individual 
ring particles is developed. Specifically, impurity levels 
and the phase of the water ice in the particles are poorly 
known (Weidenschilling , 1984), and the surface texture of 
the particles is still being debated (Kerr, 1985). 

Another feature of the model, which has received little 
attention is the collision probability, P . In the literature 

^For example, Carlson (1980) used S = 0.22 for collisions 
involving H atoms, based on experiments (Brackmann and Fite, 
1961) in which a H beam impinged upon a liquid-nitrogen- 
cooled copper plate. The beam was so dense that recombination 
of atoms to produce H 2 was a significant source of loss of 
the atoms from the beam. This experiment does not appear 
to closely approximate the physical system. 
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P c is consistently expressed as a function of the optical 
thickness, x, the definition of which is discussed in Section 
III.B. Dennefeld set P c equal to x, while recent works 
(Carlson, 1980; Ip, 1983) use 



P = 1 - e T . 

c 

This seems to be a simplification of the collision process, 

since one would expect P c to depend on such trajectory-related 

factors as velocity of the molecule relative to the ring 

particles at the location of passage, and the path length 

through the ring (Durisen (1984) states that P c will depend 

on optical thickness, and the "slant path of an ejectum 

through the rings," but does not elaborate). Also, the above 

expressions for P c do not account for multiple collisions 

0 

during one passage through the rings by a molecule. The 
existence of multiple collisions, if significantly probable, 
has implications both for density of the atmosphere (there 
are more "sticking" opportunities) , and for its spatial 
distribution . 

The aim of this paper is to elevate the theory of colli- 
sions between ring particles and atmospheric molecules to a 
level of sophistication on par with that of ballistic 



g 

Note that nowhere in Dennefeld' s model are multiple 
collisions considered. The impact of this omission is 
discussed later. 
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transport models now utilized in the study of ring atmos- 
pheres. A generalized model for determining P c is developed 
in Section II, while in Section III its applicability to 
the actual ring environment is discussed. Section IV is a 
mathematical analysis of the expression derived for P . 
Finally, in Section V the theory is applied to the special 
case of low energy, isotropic production processes (e.g., 
sublimation) to illustrate inferences that can be made on 
the spatial extent of the ring atmosphere, based on the fre- 
quency of collisions. 
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II. THEORY 



In this section, the general problem of a molecule passing 
through a translating target consisting of randomly dis- 
tributed macroscopic particles is addressed. The proba- 
bility that it collides with at least one particle is 
calculated . 

Consider a region in space of large extent, L, in two 

dimensions, and of a smaller thickness, h. This slab is 

populated by noninteracting spherical particles which are 

randomly distributed (spatially) and fixed with respect 

to the slab. Assume that the size of the spheres is limited 

to radii in the interval [r . , r << h] upon which the 

min max £ 

radii may be distributed either discretely or continuously. 

The number density of particles of size r is denoted by 

n(r).^ Finally, the slab as a whole translates in a direction 

normal to its thickness with constant velocity, v . 

P 

Now suppose that a molecule traveling along a straight 
path with speed v m passes through the slab. Assume that the 
molecule is much smaller than the smallest particle in the 
slab, and that it only interacts with the particles through 
mechanical collision. The relationship between the molecule 



7 -3 

For discrete distributions n(r) has units m . For 

continuous distributions it has units m - ^, so that n(r)dr 

is interpreted as the number of particles per m-3 in a radius 

increment {r,r+dr] . 
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and the slab can be conveniently described using direction 
cosines as shown in Figure 2(a). Notice that the x-axis is 
chosen to be parallel to v . 



z 




(a) v p ? 0 



Figure 2. Alternate 



z ' 





Frames of Reference 



Calculation of the probability of the molecule colliding 

with at least one particle as it passes through the slab, 

henceforth denoted as P , can be accomplished in the frame 

of reference shown in Figure 2(a). This procedure, however, 

is mathematically tedious (see Appendix B) . The problem is 

simplified by converting to a frame of reference in which 

v equals zero (Figure 2(b)). This approach is detailed 
P 

below . 

Suppose that the conversion to the frame of reference 
with Vp equal to zero has already been accomplished (see 
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Appendix C for details). For simplicity, let the path of 
the molecule be normal to the slab, i.e., y' equals zero, and 
n(r) be discrete (extension of the results for this case will 
be generalized to cases with continuous r and arbitrary 
molecular path) . The molecule will collide with any parti- 
cle of radius r^ which has its center anywhere in the cylin- 
drical volume (collision region) whose axis coincides with 

g 

the molecule's path as shown in Figure 3. The probability 



r . 

1 




molecular path (normal to 
the slab) 




i 

i 



Figure 3. Collision Region for Particles of Radius r 



that a collision does not occur with particles of size r^ 
is equivalent to the probability that no particles of that 



g 

Note that it was tacitly assumed that the molecule is 
essentially a point mass. This enables one to define the 
collision region as was done above. An equivalent way to 
conceptualize the collision cylinder is to let the mole- 
cule have radius r, and the particles be point masses. The 
molecule then collides with any particles it "sweeps out" 
along its path. Of course the resulting collision region 
again is a cylinder of radius r. 
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size have centers in the collision region. Denote this 
probability as P nQ (r^). Since the particles were assumed to 
not interact, the p no ( r j_) for all i are mutually independent. 
This allows one to write for discrete r: 

P c = 1 - II p no > over all r^ . (II. 1) 

i 

The problem is to find an expression for P nQ (r^). 

Under the assumptions of the model, the uniform proba- 
bility distribution governs the location of each particle, 
i.e., the probability that a particle will be found in some 
portion of the slab is proportional to the volume of that 
portion. For example, consider one particle of size r^. 

The probability that it is not in the collision region is 

L 2 h -ur 2 h 

P n°(l) - ' 

where the numerator is volume of the slab minus the collision 
region, and the denominator is the volume of the slab. Since 
the location of each particle in the slab is independent of 
the location of other particles, the probability that any 
two specific particles of size r^ are not in the collision 
region is 

L 2 h -7Tr 2 h 2 

P no (2 l * ^ > • 
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Extending this line of reasoning, it is clear that 



2 2 

L h -irr h N (r . ) 

P nn (r.) = ( 5 — ) 1 

no l _ 2, 

L h 



where : 



N (r^ ) = n (r^) L h 



is the total number of particles with radius r^ in the slab. 

P (r.) can be rewritten as 
no l 



p _ (r . ) 
no i 



nrf n (r . ) h N (r . ) 

(1 _ — i i — ) ■ L 

( N (r i ) j 



If N(r^) is large, then 



P„^(r.) = exp{-7Tr . n (r . )h} 

no i li 



(II. 2 ) 



(II. 3 ) 



Substituting into equation (II. 1 ) with N(r^) large for all 
r^ yields 



2 

P = 1 - II exp{-nr .n (r . ) h} 

c .^11 

l 

or 

2 

P c = 1 - exp{-uh Jr^n(r^)} , 



(II. 4 ) 
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recalling that r is discretely distributed, and the molecular 
path is normal. 

Extending this result to the case when the molecule's 
path makes an angle y ' with a normal from the slab is easily 
accomplished. Notice that the only difference between this 
case and that of the normal path is that the length of the 
molecular path (hence the length of the collision cylinder) 
increases from h to h/|cos y' | . Therefore, 



= 1 - exp{- 



-TTh 



cos y 



l r 2 n (r . ) } 

V 1 1 



(II. 5) 



with -r discrete and the molecular path- not necessarily normal. 

Extension to the case when r is continuous is affected 
by replacing the summation in equation (II. 5) with integra- 
tion (see the first footnote on page 9 to see the new inter- 
pretation required for n(r)). This yields 



P 

c 



_ , max 9 

1 - exp{-r— ^ , | / r n(r)dr} , 

| cos y [ 1 

min 



(II. 6) 



where r is continuous and the molecular path is not neces- 
sarily normal to the slab. 

Equations (II. 5) and (II. 6) represent solutions to the 
problem when working in a frame of reference for which 
equals zero. Often problems are not initially formulated 
in this particular frame. Appendix C shows how to convert 
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these results back to the original frame of reference (shown 
in Figure 1). Substituting from equation (C.8) yields 



1 -exp{ 



-Trh 
cos y 




(II. 7) 



for r discrete, and 



= 1 -exp{- 



-rrh 



v 



max 



cos y 



1-2 (^E.) cosa+ (^2.) / r n(r)dr} 



m 



m 



min 



(II. 8) 



for r continuous. Equations (II. 7) and (II. 8) are more con- 
venient than equations (II. 5) and (II. 6) for some applications 
It is interesting to note that while the length of the 
molecular path in the frame of reference with not equal 
to zero is h/ | cos y|» the "effective path length," pertinent 
to the collision problem, is 




(II. 9) 



which accounts for the relative motion of the molecule versus 
the ring particles. 

Equations ( II . 5 ) - ( II . 8 ) represent the fruition of this 
derivation of P . Table 1 gives a review of the tacit and 
explicit assumptions of the model. 
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TABLE 1 



ASSUMPTIONS OF THE THEORETICAL MODEL 



Assumptions 


Comments 


spherical particles 


mathematically convenient 


random spatial distribu- 
tions of the particles 


justifies using the uniform 
probability density function 


h < L 


facilitates application of 
the model to the ring-s 


r . < r < r 

min — — max 


no significance beyond 
"quantifying" particle size 


r << h 

max 


effectively excludes colli- 
sions outside the nominal 
confines of the slab 


mm molecule 


allows the molecule to be 
treated as a point particle 


molecule enters and leaves 
slab through its large 
(area of L^) faces 


facilitates application of 
the theory to the rings 


slab translates with 
velocity, v , normal to 
its thickness, h 


facilitates application of 
the theory to the rings, 


v and v constant 
p m 


mathematically convenient 


particles fixed relative 
to the slab 


mathematically convenient 


number of particles in the 
slab is large 


justifies exponential 
approximation in probability 
calculation 



Is the theoretical model utilized here relevant to condi- 
tions in Saturn's rings? This question is discussed in the 
next section. 
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III. APPLICATION TO SATURN'S RINGS 



A. APPLICABILITY OF THE THEORETICAL MODEL 
1. General Comments 



Before proceeding with analysis of collisions between 

ring particles and atmospheric molecules, the validity of the 

generalized "slab" model of Section II must be established 
\ 

vis-a-vis the actual ring environment. The essential assump- 
9 

tions of the model are: 



(1) the particles in the slab are spherical. 



(2> the radius (size) of any particle in the slab is 
in an interval [r m j_ n , r max ] where r m -j_ n is much 
greater than the size of the molecule and r max 
is much smaller than the slab thickness, h. 

(3) the spatial distribution of particles within the 

slab is random, i.e., the location of each particle 
is governed by the uniform probability density function 
with all particles mutually independent (in a sta- 
tistical sense) . 



(4) the number of particles within the slab is large 
enough for the exponential approximation (equation 
(II. 3)) to be valid. 



(5) the velocity of the slab, v , is constant. 

(6) the velocity of the molecule, v , is constant. 

(7) the slab particles are fixed relative to the motion 
of the slab as a whole. 



9 

Some of the assumptions listed in Table 1 were adopted 
to facilitate application of the model to the rings. These 
assumptions have no physical significance, so they are not 
discussed here. 
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If a "physical slab" can be defined in the rings which 
reasonably satisfies the above conditions, then the appli- 
cability of the model is verified. Note, however, that the 
physical slab is an abstract concept, since it is related to 
the passage of a specific molecule, on a specific trajectory. 

To facilitate the discussion of the above seven 
assumptions, they can be divided into three categories, those 
related to: 

(1) particle properties (assumptions (1) and (2)). 

(2) particle distribution properties (assumptions (3) 
and (4 ) ) . 

(3) velocity-related properties (assumptions (5)-(7)). 

Below each of these categories is discussed with respect to 
Saturn's A and B rings. The other rings will not be con- 
sidered, since they are not expected to contribute much to 

the ring atmosphere (owing to their low masses) . Observational 
data referenced in Section I supports this exclusion. 

2 . Particle Properties 

Individual ring particles have not yet been observed, 
thus absolute verification of the assumption of spherical 
shape cannot be accomplished. However, such a shape is 
adopted ubiquitously in the development of ring models. 
Observational data does not suggest that this traditional 
approach is invalid (Weidenschilling et al., 1984), in fact 
the Voyager radio occultation experiment (Tyler et al., 1983) 
eliminates the possibility of very elongated particles. 
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Therefore, at present the most reasonable shape to attribute 
to ring particles is spherical. 

Most ring particles have sizes in a continuous radius 
interval (Cuzzi et al . , 1984; Esposito et al., 1984; 
Weidenschilling et al . , 1984) 

1 cm < r < 5m. 

3 

Clearly, r . is much greater than the size of a H-0 or H 
molecule. The thickness, h, of the A and B rings is nominally 
not larger than 200 m, but locally may be as small as 10 m 
(Cuzzi et al., 1984). The requirement that r be much less 
than h should not create a problem for application of the 
collision theory. 

3 . Spatial Distribution of the Ring Particles 

The simplest way to describe the rings is as a 
collection of mutually colliding particles whose trajectories 
are dominated by the gravitational influence of Saturn. Under 
such conditions the notion of a completely random particle 
distribution throughout the rings is reasonable. The physi- 
cal slab could then be defined as any translating region 
containing the molecular path, provided that the velocity- 
related assumptions are valid for that specific region. 

Choose the largest such region consistent with the other 
assumptions to be the physical slab. As will be shown in 
the discussion of velocity-related requirements, this region 
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will have a considerable volume. This is desirable, for even 
if n(r) is quite small (as limited measurements have indi- 
cated (Marouf et al., 1983)) , the requirement that the number 
of particles in the slab be large is likely to be satisfied. 

In actuality the rings exhibit considerable fine 
structure (i.e., variations in optical thickness) on radial 
scales down to -1 km (Esposito et al . , 1984). Apparently 
the primary cause of this fine structure is linked to the 
gravitational effects of other bodies in the Saturnian system, 
especially the moons. Gravitational resonances are respon- 
sible for the creation of zones where stable particle orbits 
are precluded (e.g., the Cassini Division), and may generate 
density waves within the rings. The net effect of fine struc- 
ture on the collision theory is to limit the scale over which 
random spatial particle distribution is correct. This can be 
best understood by considering the following examples. 

In the first example, the slab is divided into two 
regions as shown in Figure 4. Suppose that the particles in 



t 




Note : For simplicity, 

consider the frame of 
reference with v = 0, 
and require the fnolecular 
path to be parallel to 
a slab edge. 



Figure 4. Model for Analyzing the Effect of Nonrandom 
Particle Distribution 
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each region are constrained to remain within that region; 
this represents the departure from truly random spatial dis- 
tribution of the particles. Let n^(r) equal n 2 (r) with 
r discrete. What is the probability that the molecule 
suffers at least one collision when passing through the entire 
slab? Attack this problem by reformulating it as passage 
through two smaller slabs shown in Figure 5 as shaded regions. 



Figure 5. The Two Slabs used in Reformulating the Problem 

From Section II, the probability of no collision with parti- 
cles of size r. for the slab of thickness b is 




h-b 



1 




(1 





where : 
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For all this becomes 



no (b) 



n (i 
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which for large N(r^) becomes 
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Similarly, for the slab of thickness (h-b) 



r TT (h-b) r Z , <■, 

= expt- -i — 1 tv ) r.n(r.)} 
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no (h-b) 



(III . 2) 



for n(r^)L (h-b) large. The probability of collision for the 
entire passage through the original slab is 



P c 1 P no (b) P no (h-b) 



(III. 3) 



Substitution yields 



-Tib 



P = 1 - exp{-i M1J , | J rfn (r . ) }exp{ i — — ^ j r 2 n (r . ) } , 

^ |cos y 1 f • ii ^ |cos y'I • ii 
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or 



P 



c 



1 - exp{ 



cos y ' J 



-iTh 



l r 2 n(r . ) } . 
v x 1 



(III. 4) 



This result is the same as that obtained when a random spatial 
distribution over the entire slab is assumed. Careful con- 
sideration reveals an important qualification to this conclu- 
sion. Notice that the assumption of exponential form is less 

tenable in the "nonrandom" case, since there it is required 

2 2 

that both n(r^)L b and n(r^)L (h-b) be large, while in the 

2 

completely random case n(r^)L h must be large. Extending this 
line of reasoning to cases where the original slab is further 
subdivided into regions to which the particles are confined, 
it is clear that at some point the exponential approximation 
will certainly fail. 



two regions (Figure 4), but now n^(r) is not equal to ri 2 (r) . 
Following the approach taken above, the collision probability 
is 



In the second example, the slab is again divided into 




T l r i n l (r i } }exp{ T 
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(III. 5) 
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where both n^(r^)L b and n^Cr^L (h-b) must be large enough 
to justify the exponential approximation. The collision 
theory is here valid in principle, but its application is 
more difficult. 

From this point on, the A and B rings will be assumed 
to satisfy the spatial distribution conditions. However, it 
must be kept in mind that fine structure effects may inter- 
ject inaccuracies into calculations. 

4 . Velocity-related Assumptions 

Reviewing, the three velocity-related assumptions 
are: v^ is constant; v m is constant; and the particles in the 

slab -are fixed relative to v . Each in turn is addressed 

P 

below, but the validity of the concept of slab velocity when 
applied to the ring system must first be demonstrated. As 
discussed in Section III. A. 3, choosing the physical slab to 
be as large as possible is advantageous. By means of a simple 
example, it can be shown that a single velocity well approxi- 
mates the orbital motion of all of the particles in even a 

quite large slab (this velocity corresponds to v ) . Suppose 

P 

that a slab with L equal to 10,000 km is centered in the B 
ring as shown in Figure 6. If v is set equal to the parti- 
cle velocity at C, v , how much error is inherent to the 
system? First, consider the magnitude of v . The location in 
the slab with the greatest deviation in speed relative to v c 
is either point A or B. In general the orbital speed for a 
ring particle in the Saturnian system is 
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Figure 6. A Slab with L = 10,000 km Centered in the 
B Ring 

v = /GM /R = 6159 _ [j^ns - !] . (III. 6) 

S /R 

Applying this equation at A, B, and C yields 

= 19.496 kms ^ , 

Vg = 18.587 kms ^ , 

and 



v = 19.025 kms 

c 

Thus, all of the particles in the slab orbit Saturn with a . 
speed in the range [18.587 kms ^ , 19.496 kms the greatest 
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deviation in assigning v c (equals v ) to all of the particles 
is 2.5%. Now, consider the direction of v . The point of 
greatest deviation in direction, <p , occurs at point E (the 
angle of deviation at point F is equal to that at E) . Figure 
7 gives the computation of <p ; the results given there demon- 
strate that the angular error in assigning v c to all of the 




<t> 

<P 



v T 



v 



Ex 



V 



Ey 



• ^ -l r 5000 
.. y 99 ,190 

= 5.006 xicf 2 
= 99,915 km 



■) 

rad 



=19.485 kms 
=0.975 kms 
=19.460 kms 



Figure 7. Calculation of the Maximum Angle of Deviation, -j> 



particles is quite small. The above example verifies the 
concept of a slab velocity, v , for large slabs in the B 
ring (similar results also hold for the A ring) . The concept 
applies even better for smaller slabs. 

Accepting that v^ is a meaningful concept with regard 
to physical slabs in the A and B rings, the requirement that 
it be constant can now be examined. Once defined, the slab 
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will translate as shown in Figure 8. Upon reflection, the 




Figure 8. Translation of the Slab over Time 



equals | v ^| , the direction of the velocity vector changes 
with time. The angle <p which quantifies the change in direc- 
tion is given by 



$ 




(III. 7) 



For a slab centered in the B ring (as in Figure 6) , 



<P 



8156 xlO 4 t 



(III .8) 
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where t is in seconds, <p in radians. The time required for 

-2 

a variation of 1.7453 x 10 rad (1°) is 96.1 s. This example 
suggests that for the A and B rings, the assumption of con- 
stant v is approximately correct, if the time of passage 
for a molecule through the ring is short. The requirement 
will be met for most passages (i.e., over most y), since the 
thickness of the ring is so small. To see this, then sup- 
pose that: 

(1) h = 100 m, 

(2) R = 105,000 km (molecule passes near the center of 
the B ring) , 

(3) | v | =20 kms* 1 . 10 

Over what range of y (recall Figure 2) will (j> be less than 1°; 
or phrased differently, over what range of y will the time 
of passage of the molecule be less than 96 s? In general the 
time for passage for the molecule is 



t = 



v m cos y 
m 1 



(III. 9) 



Substituting the given information, requiring t to be equal 
to 96 s, and solving for y yields 



96 



100 



20x10 | cos y | 



/ 



^The lowest kinetic energy of a stable orbit (measured 
in a frame of reference fixed with respect to Saturn) for an 
atmospheric particle which originates in the A and B rings, 
as it passes through those rings, is ~20 kms~l. 
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| cos y I = 5.21x10 , 

y = 89.997°, 90.003° . 

Therefore, in this representative case the range of y is 
0 < y < 89.997° , 90.003° < Y 1 180° ; 

virtually any y meets the requirement. 

-y 

The requirement that v^ be constant, as in the case 

-y 

of v , cannot be strictly met, since the trajectory of any 
P 

molecular orbit will have some curvature. 1 '*' However, the 
extreme thinness of the rings allows the condition to be 
approximately satisfied, so this aspect of the theoretical 
model can be applied to Saturn's rings with no problem. 

That the particles on the slab are fixed with respect 

—y 

to Vp cannot be strictly satisfied, either. Direct measure- 
ments of the relative velocities of the ring particles have 
not been accomplished, but various indirect methods yield 
consistent results: ~1 cms ^ (Weidenschilling et al . , 1984). 

Compared to the range of orbital speeds for particles in the 
A and B rings 



16.6 kms ■*" < 

(outer A ring) 



v 

P 



< 



20.3 kms , 
(inner B ring) 



Note that the gravitational interaction between the 
molecule and ring particles has been neglected, a reasonable 
assumption since the total mass of the rings is less than 10 
that of Saturn (Null et al . , 1981). 
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calculated using equation (III. 6), the relative motion of 
the particles is insignificant over reasonable passage times 
(less than or equal to ~100 s, as discussed above). Thus, 
in this regard the theoretical model suffers no important 
restriction in applicability to the rings. 

5 . Summary 

The above discussion has demonstrated that the condi- 
tions upon which the theoretical model is based, are gener- 
ally well-satisfied in the ring environment. The only 
potential source of significant concern relates to assumption 
of random particle distribution and the validity of the exponen- 
tial ...approximation (Section III. A. 3). Unlike the other assump- 
tions of the model, these cannot be unamibguously shown to 
hold for Saturn's rings. On the other hand, they cannot be 
shown to fail, and have considerable intuitive appeal. With 
this caveat noted, the theoretical model is adopted for analy- 
sis of collisions between ring particles and atmospheric 
molecules. 

B. THE UTILITY OF OPTICAL THICKNESS IN APPLICATIONS 

Having shown that the theoretical model can be viably 
applied to the atmosphere in the vicinity of the A and B 
rings, calculation of P c using the results of Section II is 
justified. Since particle radius is continuously distributed, 
and effectively confined to an interval [0.01 m,5 mj as dis- 
cussed earlier, the collision probability is (from Section 
II) 
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, r -7lh 

1 - expt-! 

^ cos y 



(II. 6) 
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c 



/ r 2 n (r ) dr } , 

0.01 



when working in a frame of reference with equal to zero, 
and is 



= 1 - exp{- 



-irh 
cos y 



' 1-2 



v 

(-E 

v m 



) cos 



V 

a + (^) 

m 



r n (r ) dr } 



0.01 



(II. 8) 



when working in a frame of reference where is not equal to 
zero. The value for h is well-constrained, while the param- 
eters v^, v , a, y t and y' depend on the location of a mole- 
cule's passage and its trajectory (both of which can be 
determined using Keplerian mechanics). If n(r) is known 
throughout the rings, then the calculation of P c using the 
above equations will be straightforward. Although the func- 
tional form of n(r) is believed to be 

n (r ) = n Q r p , p > 0 , 



where p may differ over r, its precise expression has only 
been determined for a few locations in the A and C rings 
(Marouf et al . , 1983). Fortunately, this obstacle to the 
application of equations (II. 6) and (II. 8) can be overcome 
by using the concept of "optical thickness." 
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12 

Optical thickness , x , is defined in an operational 
manner: if an electromagnetic beam traverses a layer of 

particles at an angle y to the layer normal, and the layer 
has an optical thickness of x, then the beam's intensity is 

— T / CO S U 

reduced by a factor of e . In general x will depend 

on: particle composition, size, and shape; particle density; 

layer thickness; and wavelength of the incident beam. De- 
tailed analysis of this phenomenon is beyond the scope of this 
paper (see van de Hulst (1981) and Chandrasekar (1960) for 
thorough analyses of light scattering by media consisting of 
small particles). However , the essential attributes of a com- 
plete- analysis of x, can be concisely summaried, generally 
following Cuzzi et al . (1984 ). 

For spherical particles in a layer of thickness, h, the 
optical thickness is given as 

T ( A ) = / Q e ( r , A ) rr^hn (r ) dr , (III. 10) 

where : 



A = wavelength of incident radiation. 



and 



12 

Optical thickness is often referred to as "optical 
depth." The terms are used interchangeably in literature 
concerning planetary rings. Also, since x depends on wave- 
length, one must be careful when comparing the results of 
different studies. 
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Q (r,A) = extinction cross section (a measure 

of the portion of incident energy 
removed from the beam by absorption 
and scattering) . 

When r is much greater than A , then 

1, for incoherent incident radiation. 

2, for coherent incident radiation. 

Considering visible wavelengths with respect to the predominant 
particle sizes in the A and B rings, i.e., [0.01 m,5 m] , the 
above values for Q g apply. Therefore, for the A and B 
•rings' the visible optical thickness when measured with an 
incoherent signal (e.g., stellar occultation) is 

x = / 7 Tr^hn(r)dr . (III. 11) 




Equations (II. 6) and (II. 8) can be rewritten as 

p = 1 - e ‘ T/|cos Y ' I 



(III . 12) 
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c 



1 -exp{- 



-T 



cos y 



j l-2 (^)cos 
m 



a + 



v 

(-£ 

V 

m 



r } 



(III. 13) 



Of course, if optical thickness measurements based on coherent, 
visible signals are used to compute x, the resulting expres- 
sions for P c will differ from equations (III. 12) and (III. 13) 
by a factor of 0.5 in the exponent of e. 
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Measurements of t for Saturn's rings abound. To apply 
equations (III. 12) and (III. 13) to the rings, one must only 
be cognizant of the importance of the wavelength at which 
x was determined. Keeping this in mind, analysis of P c can 
proceed . 
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IV. GENERAL ANALYSIS OF THE EXPRESSION FOR 



THE COLLISION PROBABILITY, P„ 

A. RESTRICTIONS ON THE PARAMETERS 

When studying Saturn's rings, the frame of reference 
usually adopted features Saturn as a fixed central body 
orbited by ring particles following (circular) Keplerian 
trajectories. It is also the reference frame of most general 
interest when examining the ring atmosphere, since ballistic 
transport models are based on basic orbital mechanics. For 

this reason only equation (III. 13) is analyzed in this 

. 13 

section. 

Recall that the collision probability is given as 



P 

c 



1 - exp{- 



-T 



cos 



t r 



1 1 - 2 (^£)cos 
m 



a + 



v 

m 



(III. 13) 



Clearly P c is a function of parameters related to the mole- 
cule's orbit (v m ,a,y) , and those that characterize the 

translating slab (t,v ). Note that a and y are not indepen- ■ 

P 

dent, but must satisfy 

2 2 

cos a + cos y < 1 , 



13 

That is not to say that the reference frame with slab 
particles at rest is not useful in some practical applications 
(as is shown in Section V). Rather, although the collision 
problem may be solved in either frame of reference, the prob- 
lem of calculating molecular trajectories is ultimately ap- 
proached in the frame with v not equal to zero, requiring 
conversion to that frame. ^ 
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since they are direction cosines. An alternative coordinate 
system which removes this interdependence specifies a mole- 
cule's path in terms of y and 9, henceforth referred to as 
the "orbital inclination" and "orbital azimuth." 9 is here 
defined as the angle between the projection of v m in the 
x-y plane (the ring plane), and the x-axis (i.e., the direction 
in which the slab translates) . Figure 9 shows the relation- 
ship between 9 and the direction cosines, the mathematical 
expression of which is given in Appendix D. Employing the 




Figure 9. Definition of the "Orbital Azimuth," 9 

results of Appendix D, the collision probability can be 
rewritten as 



-T 



P c = l-expt-j co3 




V 



V 



1-2 (— £) sin y cos 9 + (-2-) } 



v 



m 



v 



(IV. i: 



m 



Equation (IV. 1) proves more useful than equation (III. 13) 
in the subsequent analysis. 
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Before examining the expression for P c in detail, physi- 
cal considerations which limit the range of values for the 
parameters are outlined below. Following this, attention can 
then be turned to parameter values of most practical interest 
in the study of the atmosphere of Saturn's rings. 

The optical thickness, x, of Saturn's rings has been 
studied extensively. Recall from Section III.B that the 
measured value of x is dependent on several variables, most 
importantly wavelength and coherence of the incident radiation 
and the nominal size of the ring particles at the location 
where the electromagnetic beam is incident upon the ring plane 
As shown in Section III.B, measurements based on incoherent, 
short wavelength (e.g., visible) radiation yield values of 
x that can be most easily used in collision theory (see 
equations (III. 11) and (III.13)). A stellar occultation 
experiment onboard Voyager 2 provided results that meet the 
above requirements (Esposito et al., 1984). Optical thick- 
ness was found to fluctuate radially, varying from: 

(1) -0.1 to ~0.4 in the C ring. 

(2) -0.7 to -2.5 in the B ring. 

(3) -0.1 to -0.2 in the Cassini Division. 

(4) ~0.4 to ~1.0 in the A ring. 

Average x over the A-C rings is given in Appendix A. Sum- 
marizing this information, the nominal optical thickness is: 

(1) 0.1 for the C ring and the Cassini Division. 

(2) 1.5 for the B ring. 

(3) 0.5 for the A ring. 
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Subsequently, in further analysis of Saturn's ring atmosphere, 

T will be limited to the interval 

0.1 £ x <_ 2.0. 

The slab speed (which is essentially the ring particle 

speed in the vicinity of the molecule's path through the 

ring) and the molecular speed appear as a ratio, v^/v^, in 

the expression for P . Call v /v the "speed ratio." Esti- 

c p m ^ 

mates of f^O molecule ejection speeds, relative to the source 
particle, for several production processes are: 

(1) ~0.3 kms ^ for sublimation (Dennefeld, 1974). 

(2> ~1 kms for meteroid bombardment (Dennefeld, 1974). 

Estimates for H atom ejection speed are: 

(1) < 3 kms 1 for meteroid bombardment (Morfill et al., 

1983) . 

(2) < 10 kms ^ for photodissociation of I^O (Durisen, 1984). 
Based on the above values, estimates for the speed ratio 
range for molecules newly ejected from particles in the A 

and B rings are 



v 

0.95 < - 2 . < 1.06 

— v — 

m 



for H^O molecules, and 



v 

0.40 < -2 < 1.49 

— v — 
m 
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for H atoms, noting that is the sum of and the ejection 
velocity, v e j • The eccentricity of the subsequent orbit, 
and the effect of collisions complicate precise specification 
of general limits on the speed ratio, but the above should 
suffice as guidelines. Therefore, further analysis of P c in 
this paper will be limited to a range 



v 

< -E < 

- v m - 



2 . 



By definition, y and 9 are restricted to [0,7r] and [ — tt , tt ] , 
respectively. Considering equation (IV. 1), it can be seen 
that - 



P ( 9 ) = P (-6 ) , 

c c 



and 



P c (y ) = P c ( tt — y ) «--»■ P c (jTT-y) = P c (j7t + Y) • 

Noting these symmetries, no information is lost by limiting 

subseauent analysis of P to 

J c 

0 < 9 < it , 



and 



0 < y 



< 
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Though P c will be analyzed over the above intervals, some 
physical restrictions on values for orbital inclination and 
azimuth are discussed in the following paragraphs. 

When a given production process is considered, the initial 
limits on y can easily be determined as shown in Figure 10. 

For example, the y-range for ^0 is: 



z z 




Figure 10. (a) General Relationship Between v e j , v m , 

and Vp. 

(b) Case for Minimum Value of y 



(1) 89° 

( 2 ) 86.6 

Similarly, 

(1) 79.6 

53.0 



< y $ 91° for 

° < y < 93.4° 
for H atoms , y is 
° < y < 100.4° 



° 6 y $ 



sublimation . 
for meteroid impact, 
initially confined to: 



( 2 ) 



127.0° 
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for meteroid impact. 

for photodissociation of 1^0 
on a particle surface. 



As collisions occur between molecules and ring particles, 
the range on y will change. If a collision is assumed to be 
inelastic, then upon impact the molecule will rebound with 
reduced kinetic energy relative to the particle. The overall 
effect is that y should go to as more collisions occur. 

For atmospheric molecules originating on a ring particle, 
the initial orbital azimuth will certainly be "forward- 
directed , " i . e . , 




< 6 



< 




/ 



as long as the ejection speed is less than the particle speed. 

— >■ 

This again follows from the fact that v is the vector sum of 
J m 

—>■ — >■ 

v and v .. Subsequent collisions should not alter this 
P S3 M 

assertion. 



B. FUNCTIONAL DEPENDENCE ON OPTICAL THICKNESS, T 

P c monotonically increases with increasing t, as can be 
clearly observed in equation (IV. 1). The rate of increase 
is given by (using equation (IV. 1)) 




i r \ / 1-2 (-E-) sin y cos 

cos y V v 1 

i mv m 



v 

e + (^) 

m 



exp{- 



-T 



cos y | 




Here the effect of the other parameters on the rate with 
which P c exponentially approaches one is explicitly displayed. 
Appendix E contains examples of the collision probability's 
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dependence on x given selected values for the speed ratio, 
orbital ' inclination , and orbital azimuth. 

The physical interpretation of this functional behavior 
is simple. Recall that for measurements at appropriate 
wavelengths , 



T 



2 

/ xr n(r)hdr 



If h is held constant, then higher values of x imply greater 
number density, n(r). Intuitively, one would expect the 
collision probability to increase with increasing number 
density (all other things held constant) , which is substan- 
tiated by the mathematical expression for P c . Similarly, 
if n(r) is held constant, increasing x implies that h must 
increase, the length of a molecule's path must increase, 
thus increasing its exposure to collision. 



v 

C. FUNCTIONAL DEPENDENCE ON THE SPEED RATIO, -£ 

v 

m 

The behavior of the collision probability as a function 
of speed ratio is very interesting. Figure 11 gives a repre- 
sentative example of the variation of P c as a function of 
Vp/v m (see Appendix E for the results of more example calcu- 
lations. For 9 less than iir , a minimum for P occurs at 

2 c 

Vp/v m greater than zero. Why does this happen? 

Some insight into this phenomenon can be gained by 



considering the partial derivative of P c with respect to 

v /v , i . e . , 
p m 
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Figure 11, 



Collision Probability as a Function of the 
Speed Ratio for t = 1.0 and y = 1/4 it and 
Various Values of 9 
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With all variables except v /v held constant, P has a 

p m c 

minimum at 

v 

^ = sin y cos 0 . (IV. 3) 

m 

From Appendix D, the right hand side of equation (IV. 3) is 
identically equal to cos a, so expressing that in terms of 
the direction cosine gives 

Vp = v m cos a ; (IV. 4) 

a minimum of P c occurs when the component of in the direc- 
tion of v p is equal in magnitude to | v p | . Notice that relative 

to variation in v /v , P cannot assume a minimum value for 

p m c 

speed ratio values greater than one. 

Since the result of the above mathematical analysis may . 
be difficult to reconcile with intuition, an example is 
given to clarify its physical meaning. Suppose that a blind 
mouse lives in a hole in the middle of a single-lane, one- 
way road. Randomly spaced steamrollers traverse this road, 
all moving with speed v . The mouse, smelling cheese in 
another hole a few meters down the road, dashes for that 

hole with speed v (in the same direction as the steamrollers). 
L m 

Given v , for what v will the mouse have the best chance 
m p 

to reach the cheese unharmed, neither being overtaken by a 
steam roller, nor running into one? The answer is that the 
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probability of collision is least (in this case zero) when 



v equals v . For v greater than v , steamrollers will 
tend to overtake the mouse; as gets very large the proba- 
bility of collision goes to one. For v^ less than v m , the 
mouse will tend to run into a steamroller. The mouse's 
situation is analogous to that of a molecule passing through 
a ring. The probability of collision for the molecule is 
minimized when it "goes with the flow" of the particles as 

much as is possible. This minimum value of P will not in 

c 

general be equal to zero for molecules passing through the 
ring, as it is for the mouse. This is the case, because 
the molecule crosses the ring plane of some angle (speci- 
fied by y ) • 

A final comment on the influence of v /v on P is that 

p m c 



as 



v 

vm 

regardless of the values of the other variables. For v^ 
very large, the ring loses its discrete character becoming 
essentially a continuum* collision becomes inevitable. 

D. FUNCTIONAL DEPENDENCE ON ORBITAL INCLINATION, y, AND 
ORBITAL AZIMUTH, 6 

y and 9 specify the path of a molecule through the 
ring. Figures 12 and 13 show two examples of the 
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Figure 13. P c as a Function of y and 0 for t = 1, ^ = 2 
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functional dependence of P c on orbital inclination and 
azimuth. 

First, consider the behavior of P c with respect to 9. 

The partial derivative with respect to this variable is 

3P v v 

~W = T 6os v I ( V } sinYsine [1 ~ 2 ( v } sin Tcose 
1 1 1 m m 

v ~l/2 / v v 0 

+ (-£■) ] exp{-i — r \ / 1~2 (-^) sinYcos0+ (-£) } . 

v cos y V v ' v 

m 1 1 ' m m 

(IV. 5) 

From equation (IV. 5) it can be deduced that P c is minimized 
when - 9 equals zero and maximized when 9 equals it (when all 
other variables are held constant). As expected, both Figures 
12 and 13 clearly confirm this result. Also, notice how 
rapidly P c approaches its maximum value with increasing 0. 

The story of the blind mouse from Section IV. C provides the 
physical interpretation for this behavior: all other param- 

eters held constant, collision probability is least when 
going with the flow of particles (or steamrollers), and is 
greatest when moving against them. 

The variation of P^ with y is a little more complicated. 
The partial derivative of P c with respect to y is 

3P c v v 2 i/2 

-x — = [t secY tan y (1-2 (-^) sin ycos 9 + (— ) z ) 

dY V V 

m m 

V V v ? “1/2 _ I V V ~Z 

- t (-E-) cos9 (1-2 (-E-) sinYCOs6-K-^-) A ) ] exp{n — r \ / 1-2 (-£) sinYcosG-K-^-) } 

m m m 1 1 1 V m m 

(IV. 6) 
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From equation (IV. 6), the condition for minimum P , when 
all variables except y are held constant, is 



sec y tan y 
2 

1 +2 tan y 



v 

(— £) cos 

v m 

v 

</> 

m 



1 + 



6 



2 



(IV. 7) 



This imposing expression has a simple physical meaning: 

(1) if v_/v m is less than one, P c is minimized when 
the_magnitude of the component of v m in the direction 

of v is equal to v . 

P M P 

(2) if Vp/v m is greater than one, P c is minimized when 
the^magnitude of the component of v„ in the direction 
of v m is equal to v m . 

Once again, the blind mouse analogy may be helpful in under- 
standing this behavior. 

Before leaving this section, a few words are in order 
to clarify possible concern over the point of "undefined" 
collision probability shown in Figure 12. This anomaly only 
occurs when v p/ v m is equal to one. For this particular speed 
ratio, when 0 is zero, and y is tt/2, then the molecule is 
moving in concert with the ring particles. Thus, if it is 
not initially in contact with a particle, it never will be. 

On the other hand, if initially in contact with a particle, 
it will never separate from that particle. Therefore, the 
notion of a collision probability is meaningless under this 
unique circumstance. 
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V. DISCUSSION 



A. GENERAL COMMENTS 

To effectively exploit the collision theory developed in 
this paper, one generally must couple it with a ballistic 
transport model which provides values for the speed ratio, 
orbital inclination, and orbital azimuth. More accurate 
estimation of atmospheric density will result, although 
improvements in this area are limited by the lack of knowledge 
of details of a molecule-particle collision (i.e., the stick- 
ing coefficient described in Section I) . The frequency of 
collisions also affects the spatial extent of the ring 
atmosphere. The trajectory of a molecule is initially 
determined by the nature of the production process, while its 
evolution is dictated by collisions. 

Consideration of the special case of low energy, iso- 
tropic production processes (e.g., sublimation) allows one 
to apply the collision theory in making inferences on the 
spatial extent of the ring atmosphere. 

B. APPLICATION TO AN ISOTROPIC, LOW ENERGY PRODUCTION 

PROCESS 

Define an isotropic, low energy production process to 
be a mechanism which ejects molecules from ring particles 
such that: 

(1) the direction of ejection relative to the particle 
is random. 
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(2) v . << v , which implies v /v = 1. 
ej p ^ p' - m 

The orbit of a molecule, prior to collisions with ring 
particles, will be very nearly circular with a radius equal 
to that of its parent ring particle. The orbital plane 
will be only slightly inclined with respect to the ring 
plane, e.g., for sublimation the orbital plane is inclined 
at most 1° (see Section IV. A) . 

When a collision does occur, two outcomes are possible: 
the molecule "sticks," thus removing it from the atmosphere; 
or the molecule rebounds, assuming a new trajectory. In 
the latter case, what can be said about the new orbit? 

Since, the collision will surely be inelastic to some degree, 
the speed of the rebounding particle relative to the ring 
particle is less than its incident speed. The speed ratio 
remains almost equal to one, implying a circular orbit. 
Following the approach used in Section IV. A, especially 
Figure 10, the primary effect of the inelastic collision is 
to reduce the inclination of the orbital plane. Thus, 
collisions act to reduce the spatial extent of the ring 
atmosphere; the details of this reduction being determined 
by the frequency of collisions, and the actual inelasticity 
of an individual collision. 

The more frequent collisions are, the more closely con- 
fined to the vicinity of the ring will be the ring atmosphere. 
The collision theory can be used to examine collision fre- 
quency. Specifically, the average probability, as a function 
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of time, that a molecule has not suffered at least one 



collision is calculated below. 

When a molecule is initially ejected, it may collide 
with a ring particle before ever leaving the ring plane. 

To calculate the probability that it collides with at least 
one ring particle, adopt a frame of reference in which 
the particles are fixed. The probability that the molecule 
is ejected from a particle in dz' is (see Figure 14) 

P(z')dz' = pj- dz ' , 0 £ z ' <_ h , ( V . 1 ) 

while the probability that it is ejected into solid angle 
dft is 

P(Q)d£3 = j siny' dy 1 , 0 <_ y 1 <_ tt , (V.2) 

z ' 



l 




Figure 14. Ejection of a Molecule from a Particle at z' 
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since the production process is isotropic. The probability 
density for a molecule ejected at y' from z' colliding with 
a ring particle before leaving the ring 



! 





( V. 3 ) 




which follows directly from equation (III. 12), when the 
effective slab thickness is taken into account. The aver- 
age probability of at least one collision for a molecule 
produced by an isotropic process is found by integrating over 
the product of equations (V.l), (V.2), and (V.3): 



Making a change of variables, x = cos y' and x = - cosy' 
in the first and second integrals, respectively, then inte- 
grating, yields 




P 



/ / ihsin y' [l-exp{ 

0 0 



-t ( h-z ' ) 



t-] dz 'dy ' 



c 



h cos y 




P 



c 



! 1 

i - k I X(1 
T 0 




( V. 4 ) 
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Therefore, the average probability that a molecule initially 
(i.e., t = 0) escapes the ring without collision is 



P esc = 1 " p c = k /q X(1 " e T/X ) dx ( v - 5 ) 

Table 2 gives P for values of x of interest for the rings. 

Notice that there is a significant chance of collision immedi- 
ately after a molecule is ejected; suggesting that the spatial 

extent of the atmosphere is biased toward being closely 
confined to the ring even before subsequent evolution of the 
system is considered. 



TABLE 2 

AVERAGE PROBABILITY THAT A MOLECULE INITIALLY 
(t = 0) ESCAPES THE RING PLANE WITHOUT 
COLLIDING WITH RING PARTICLES 



T 


P 

esc 


i — 1 
• 

o 


0 . 837 


0.5 


0 . 556 


o 

• 

r— 1 


0.390 


1.5 


0.296 


2. 0 


0.235 



During each orbital period, a molecule in a circular 
orbit traverses the ring two times. If that molecule has 
not collided with a ring particle, y' retains its initial 
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value for each subsequent passage of the rings, and the 
probability that it does not collide with a ring particle 
during each period is exp{ 



_2t - 14 

" • The probability that 



i cos Y i 

the molecule is ejected into solid angle dft without colliding 



with a ring particle is: 



c 2£il(i-exp{ - , o<y' <y- 

x ^ cosy — 2 



isiny'[/ c-(l-P (y 1 ,z') Jdz'jdy' 
0 



^siny '■ 






If is the number of completed orbital periods, then the 
probability that the molecule has not collided with any 
particle as a function of time, given in terms of T , is 

ir 



P (T ) 
esc p 



r 1 • - 

2 s my 



— ( 1-exp { — f} , 0 <y' < i-TT - 

x ^ cosy - ' 2 



iy' <-* 



_2tT d 

>exp ( "coiv ' ' 



or 



P (T ) 
esc p 



1 ( 

= — xe 

T 0 



- 2 iT /x 



p ' Yl-e‘ T/x )dx 



(V. 6) 



14 

This result can be easily understood by considering 
a slab of thickness 2h. Extending this reasoning, if T D 
periods have been completed, the probability of no collision 
for a molecule initially escaping is exp{-2xT / | cosy ' | . 
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Table 3 gives the probability that a molecule has remained 
free of collisions , as a function of time (given in orbital 
periods) for various optical thicknesses. These results 

TABLE 3 

VARIATION OF P WITH T AND T 
esc p 



P 

esc 



T 

P 


T = 0.1 


T = 0.5 


T = 1.0 


T = 1.5 


t = 2.0 


0 


0.837 


0.556 


0.390 


0.296 


0.235 


1 


0.519 


0.106 


0.021 


0.005 

, -4 


0.001 


2 


0.357 


0.028 


0.002 

-4 


1.56x10 


-o 


3 


0.255 


0.008 


1.94x10 


~0 


~0 


4 


0.186 


0.002 


~0 


~0 


~0 


5 


0.138 


0.001 


-0 


~ 0 


~0 



strongly imply collisions are quite frequent for isotropic, 
low energy production processes. This in turn implies that 
the ring atmosphere associated with this process will tend 
to be more concentrated in the vicinity of the ring than 
analysis of the distribution of particles prior to collision 
would indicate. 

C. ON MULTIPLE COLLISIONS 

The collision probability theory developed in this thesis 
cannot be used to precisely calculate the probability of 



63 



multiple collisions, i.e., the chance that a molecule col- 
lides with two, three, etc., particles during a single 
passage through the rings. Computer-simulation models may 
provide the best way to approach this problem. 

Based on the results of the preceding section, the 
probability of at least one collision is quite high in the 
A and B rings. This suggests that multiple collisions could 
be a significant factor in analysis of the ring atmosphere. 
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VI. CONCLUSIONS 



In this paper, a model for collisions between ring 
particles and molecules in the ring atmosphere has been 
developed. The resulting expressions for the probability 
of at least one collision when a molecule passes through 
the ring are 



P 

c 



1 - exp {- 



- T 



cos y 



(III. 12) 



for a- frame of reference where the ring particles are fixed, 
and 



P = l-exp{- 



-t 

cos y I 



v 



1-2 (-£) sinY cos 9 + (— ) } 



v 



m 



v 

m 



(IV.l) 



for a frame of reference fixed on Saturn in which the parti- 
cles are seen to orbit the planet. They are more complete 
than those currently employed in analysis of the ring atmos- 
phere, since specifics on molecular trajectory and the rela- 
tive velocity of a molecule with respect to the ring parti- 
cles are included. Thus, this analytic collision theory is 
compatible with computer-simulated ballistic transport 
models, since they include similar effects. 

The frequency of collisions influences both the total 
population of molecules in the ring atmosphere, and their 
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spatial distribution, as shown by Dennefeld (1974). Since 
knowledge of molecular sources is still evolving, and under- 
standing of the dynamics of an individual collision (e.g., 
Dennefeld' s "sticking coefficient") is in its infancy, 
improved estimates of molecular populations cannot be made 
based on this revision of collision theory, alone. Infer- 
ences on spatial distribution can, however, be made. When 
a low energy, isotropic production process is considered for 
optical thicknesses characteristic of the A and B rings of 
Saturn, collisions are found to be quite frequent. This 
implies an atmosphere closely confined to the vicinity of the 
rings. Extending this result to more energetic production 
processes (e.g., see Section III .A) implies that the ring 
atmosphere will be toroidal in shape. Thus, assuming a 
toroidal shape will probably be more productive than assum- 
ing spherical shape in future studies of the ring atmosphere. 

As a final note, the collision theory used here in con- 
nection with the ring atmosphere should be useful in another 
related application. Ballistic transport of macroscopic 
particles produced during meteroid impacts may account for 
some of the radial structure in Saturn's rings (Ip, 1983; 
Durisen, 1984). These particles form a "chip halo" in the 
vicinity of the rings. As long as the chips are small rela- 
tive to ring particles, the model developed here applies 
for calculation of collision probabilities. 
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APPENDIX A 



RING NOMENCLATURE AND DIMENSIONS 

Saturn's ring system consists of seven distinct rings. 

The A, B, C, D, and F rings are characterized by distinct 

boundaries. Their nominal thicknesses are not more than 

-200 m (Cuzzi et al . , 1984), and may well be less than 5-10 

m (Bridges et al . , 1984). The E and G rings are diffuse 

structures with poorly defined boundaries and nominal thick- 
3 2 

nesses of ~10 and 10 km, respectively (Cuzzi et al . , 

1984 \. Table Al is a summary of the radial structure of the 
ring system. Note that the F and G rings are quite narrow. 

TABLE Al 

RADIAL STRUCTURE OF SATURN'S RINGS 



Boundaries 



Ring/Region 


R 

s 


km 




mass 




D 


-1. 11-1. 235 


66,970-74,510 




? 




C 


1.235-1.525 


74,510-92,000 


2 


-9 

x 10 


Ms 


B 


1.525-1.949 


92,000-117,580 


5 


00 

1 

o 

i — ! 

X 


Ms 


Cassini 

Division 


1.949-2.025 


117 , 580-122,170 


1 


X 10 9 


Ms 


A 


2.025-2.267 


122,170-136,780 


1 


00 

1 

o 

p — 1 

X 


Ms 


F 


2.324 


140,180 




7 




G 


2.82 


170,100 


1 


. rt -17 
X 10 


Ms 


E 


00 

1 

00 


181,000-438,000 




7 





2 6 

R = 60,330 km (Saturn's radius) M = 5.685x10 kg (Saturn's 

s mass) 
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The Cassini Division is a well-defined gap between the A and 
B rings, which is by no means "empty." 

The optical depth, t, is a dimensionless parameter 
directly related to mass density. Table A2 gives average 
optical thickness for the main rings. Optical thickness in 
the other rings is much lower. Since optical thickness 
increases with increasing mass density, it follows that the 
A and B rings should be the most significant sources of a 
ring atmosphere. 



TABLE A 2 

AVERAGE OPTICAL THICKNESS OF THE MAIN RINGS BASED 
ON STELLAR OCCULTATION OBSERVATIONS FROM 
VOYAGER 2 (ESPOSITO ET AL . , 1984) 



REGION 


BOUNDARIES (R ) 
s— 


~AVE 


inner C 


1.24-1.39 


0.08 


outer C 


1.39-1.52 


0.15 


inner B 


1.52-1.66 


1.21 


middle B 


1.66-1.72 


1.76 


outer B 


1.72-1.95 


1.84 


Cassini Division 


1.95-2.02 


0.12 


inner A 


2.02-2.16 


0.70 


outer A 


2.16-2.27 


0 . 57 



R = 60,330 km 

s 
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APPENDIX B 



DERIVATION OF P WHEN v IS NOT EQUAL TO ZERO 
c p 



In Section II of the main body of this paper, an expres- 
sion for the probability that a molecule passing through a 
ring collides with at least one ring particle, P c , was 
derived in a frame of reference in which the ring particles 
are fixed (i.e., equal to zero). Here it will be demon- 
strated that the same expression for P c results when the 
derivation is made in a frame of reference where the slab 

-y 

translates with v^ not equal to zero. 

The model upon which this derivation is based is the 
same as that described in Section II. Reviewing, a molecule 
passes through a region in space which contains randomly 
distributed (spatially) spherical particles. This region, 
the "slab," is of large extent, L, in two dimensions- and of 
relatively small thickness, h. The entire slab moves through 

-V 

space with constant velocity, v , in a direction normal to 

—V 

the slab's thickness (v is measured relative to some iner- 

P 

tial frame of reference). The radius, or "size," of each 

particle is restricted to an interval (r . ,r ] with r 
c mm max mm 

much greater than the nominal radius of the molecule, and 

with r much smaller than h. The size distribution of the 
max 

particles may be either discrete or continuous. The number 
of particles in the slab is further assumed to be large. 
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Finally, the number density as a function of size is denoted 
by n(r), defined as the number of particles: 

(1) of radius r per unit volume (discrete size 
distributions) . 

(2) in the radius interval [r,r+dr] per unit volume 
(continuous size distributions) . 

What is the probability that the molecule collides with 

a particle in the slab? In Section II this question is 

answered with relative ease by converting to a frame of 

reference in which the particles are at rest. Below, no 

such simplification is made. 

— >- 

When v is not equal to zero, the notion of cylindrical 
P 

"collision regions," used in Section II, loses its validity. 

Another approach must be taken. Suppose that the slab is 

subdivided into regions of very small thickness. Ah, and 
2 

finite area L (see Figure Bl) . If the probability that the 









'///'/”//// / 7 7/ /-'W 


h 

jL 





k — l 



Z Ah. 



Figure Bl . Edgeview of the Slab and a Representative 
Subregion 
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molecule does not collide with particles of size r whose 



centers lie in a specific sublayer can be calculated (call 
this P nQ (r,Ah^)), then the total probability of no colli- 
sion with particles of size r will be 



P no (r) = " P no (r ' Ah i ) 

1 



(B.l) 



since all of the small sublayers are mutually independent . ^ 
The collision probability, following Section II, is 



P = 1 -n P (r.) = 1 - n H P (r.,Ah.) ( B . 2 ) 

c no i . . no i i 

D J D i J 

for discretely distributed r. 

The problem now is to find an expression for P nQ (rj,Ah^) . 
First, the case of discrete r is solved from which the con- 
tinuous solution can be deduced. Adopt a point of view in 
which the molecule has radius r, and the particles are point 
masses (on reflection this is analogous to the real situation 
as far as collision probabilities are concerned). Then, 

P (r.,Ah.) can clearly be seen to be equivalent to the 
probability that the volume in the subslab "swept out" by 



15 

Two observations on equation (B.l) are in order. First, 
there are infinitely many factors in the product over i. 
Second, P no (r,Ahj_) must be the same for all i, since in this 
model the statistical properties of each subregion are iden- 
tical. Equation (B.l) may be difficult to interpret at 
this point, but its meaning will be clarified later in this 
appendix . 
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the passing molecule has no particles centered within it. 

The probability of no collisions with particles of size 

r. in subslab Ah. is: 

3 i 



p no <r j' 4h i> 



L 2 Ah.-(Vol. Swept Out) n(r.)L 2 Ah. 
[ 1 i 3 2 



L 2 Ah. 

l 



(B . 3 ) 



since the particles are assumed to be randomly distributed 
throughout the slab (see Section II for a more detailed 
explanation of this assertion). In general, the mathematics 
involved in finding this "volume swept out" is extremely 
tedious. To demonstrate the method, the special case of y 
equal to zero or tt (molecular path normal to the plane of 
the slab) is solved below. 

On what portion of its path can a molecule potentially 
collide with a ring particle of size r_. in Ah^? Consideration 
of Figure B2 shows that the pertinent part of the path is of 
length 2r ^ . If time, t, is set equal to zero when the mole- 
cule initially encounters the subslab (Figure B2 (a) ) , then 
at any t such that 



2r . 



0 < t < 

— — v 



m 



the molecule will collide with any particle is a disk of 
radius (call this the "interaction radius") 



R I (t) 



2 2 X / 2 
(2rv m fc - V > 



(B . 4 ) 
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Molecular 

Path 




Figure B2. Relationship Between the Molecule and a 

Subslab When Collision with Ring Particles 

Centered in the Subslab is: (a) First 

Possible (t = 0); (b) At Any Time, t, 

Such That 0 < t < 2r./v 

D m 



where R^ is defined as shown in Figure B2(b). As the inter- 
action radius evolves, the particles in the subslab move 
with speed v . Considering these two factors together, it 
can be deduced that an oblong "interaction region" in the 
subslab must be particle-free, if no collisions are to occur. 
To actually calculate this volume, first relate R^ to the 
distance that the slab moves in time t, call this ; clearly 

D ( t ) = v t . (B.5) 

I P 



Substituting from equation (B.5) into equation (B.4) yields 



VV 



v v _ 1/2 

(2r (^i)D I - (^) 2 Dj) 



(B. 6) 
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Figure B3 is a sketch of equation (B.6) . The oblong area 
in Figure B3 (bounded by the graph for and its reflection 
across the D^-axis, which is shown by a dashed line) cannot 




Figure B3. as a Function of 

represent the projection of the interaction region into the 
plane of the subslab. Consideration of instantaneous inter- 
action region, represented by the circle in Figure B3 , illus- 
trates that the actual interaction region is larger. Knowing 

Rl(Dj), the function enveloping the interaction region can 

* 

be found. Call this function R^ . 

Consider Figure B4 ; it shows one-quarter of the inter- 
action volume's projection into the plane of Ah^. At a given 
time t', D' equals v t', and the interaction radius is R ' . 

The value of R" (D ) for any is 

M [R I (D ]; ) [ 2 -[D^-D I ] 2 ) 1/2] ; lDj-D'1 <R^(D') . 

R^(Dp = | (B . 7 ) 

/ 0 ; otherwise . 




R*(D^) is the maximum value of R^!(D.j.). Setting the first 
derivative of R^ with respect to equal to zero, and 
solving for gives 



D i 



D + r . (— ) 

1 3 V p 

[1 + (^) 2 1 

P 



(B . 8 ) 



when equation (B.6) is substituted into equation (B.7). 
Substituting from equation (B.8) into equation (B.7) gives 
the maximum value of R^(D^) for a given D , i.e.. 



R I* (D I } 



[(>{ 

P 



2 , V m, , 0 _ _2 , V m N 

r j ( — 1 1 +2r j D I -°I ( — 1 1 1/2 

d_£ E_>] 

v o 

1 + (— ) 2 
v 

P 



(B . 9 ) 
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Using equation (B.9), the region in the sublayer that must 
be devoid of particles of size r^ if no collision is to 
occur between the molecule and such particles, has 



Volume 



(Ah i ) 4 / 



V 

r, (rr^) 
m 



3 v. 



R*(D ) dD 

D 1 1 

I (min) 



(B . 10 ) 



where 



I (min) 



v v - 1/2 

((-£) - {1 + (-£) } ] r . 

v m v m ^ 



(B.ll) 



When the integral in equation (B.10) is evaluated. 



9 v 0 1/2 

Volume = rrr.{l + (^) } Ah . , 

3 m 1 



(B . 12 ) 



and substituting into equation (B.3) gives 



? 7 V 7 l/ 2 7 

L z Ah.-^r {!+(-£) } Ah. n(r.)L Ah. 



P no (r j' Ah i ) ' >■ 



i ] v 



m 



L 2 Ah . 

l 



l 3 
-] 



(B. 13a) 



Equation (B.13) can be rewritten as 



p no< r :- ah i ) = 11 - 



? v _ o 1/2 

irr ,n(r . ) {1 + (^) } n (rjL 2 Ah, 



m 



-] 



L n(r.) 



(B . 13b) 
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Summing over i yields the probability that no collision 
occurs in the entire salb , i.e., 



P (r . ) 
no 3 



[1 - 



2 v d 2 1/2 

7 Tr j n(r j ){1 +( % } } n(r.)L 2 h 

L 2 n(r_.) 



(B . 14a ) 



Multiplying the second term in the brackets in equation 
(B.14)- by h/h yields 



P (r . ) 
no 3 



= [1 - 



, v - 1/2 

*r 2 n(r >{1 + (UL, } h n(x ^ 

m ] 3 

L 2 n (r ) h 



(B . 14b) 



but the total number of particles of size r ^ , n(r^)L h, is 
assumed to be large, so equation (B.14b) can be rewritten 
in exponential form as 



P (r . ) 
no 3 



, v 0 1/2 

exp{-7Tr^n (r . ) {1 + (-2-) } h} 

- 1 - 1 v m 



(B . 14c ) 



Therefore, the collision probability is (from equation (B.2)) 



P 



c 



1 



n 

j 




2 1/2 
2 } h} 



/ 



p 

c 



V 



1 - exp{-TT (l r n (r . ) ) {1 + (^) 
j - 1 ^ v m 




h} , 



(B.15a) 



for r. discrete, or by extension 
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(B . 1 5b ) 



? • v 9 1/2 

P„ = 1 - exp{-7T (/ r n(r)dr) {1 + (-£) z } h} , 

w V 

m 



for continuous r. 

Comparison of equations (B.15a) and (B.15b) with equations 
(II. 7) and (II. 8) (with y equal to zero), reveals the equiva- 
lence of solving the collision probability problem in two 
different f rames-of-ref erence . That the results are the same 
is not surprising. What is interesting is the relative 
difficulty of the two approaches. 



78 



APPENDIX C 



FRAME OF REFERENCE CONVERSION 

In Section II of this thesis, the collision probability, 

P , is calculated in a frame of reference in which the ring 

particle velocity, v , is equal to zero. When coupling the 

P 

collision theory with ballistic transport models, a more 
convenient frame of reference is a localized cartesian sys- 
tem fixed in the vicinity of the intersection of the mole- 
cule's path and the ring. The x-axis is in the direction of 
motion of the ring particles; the y-axis is directed radially 
outward from Saturn; and the z-axis is normal to the ring 
plance (see Figure Cl). Here the direction cosines of a, 6, 



z 




Figure Cl. A Coordinate System Fixed at a Specific 
Location in the Rings 
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and y quite accurately describe the path of the molecule 
through the ring which is essentially linear owing to the 
great radius of curvature of the molecular trajectory 
(~10~* km) relative to the extreme thinness (—10 ^ km) of the 
rings . 

Let the frame of reference, for which the ring particles 
in the region through which the molecule passes are at rest, 
be called the "primed" system. The simplest coordinate 
system that can be used to describe the molecule's path in 
this frame of reference will have the x'-, y'-, and z'-axes 
in the same direction as the corresponding axes in the un- 
primed system with the y ' -axis in the direction of v^ (see 
Figure C2) . a', 8', and y' yield direction cosines describing 

the molecules path in the primed system. 

z z ' 





Figure C2 . Direction Cosines in the Two Coordinate Systems 



In Section II the collision problem is attacked using 
the primed system. The expressions for P c derived there 
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(equations (II. 5) and (II. 6)) contain a factor of r cos y'J ^ 
An expression for this factor in terms of a, 3, y, v^, and 
Vp must be found to represent P c in the unprimed system. 
First, note that: 





(C.l) 



gives the relationship between the molecular velocities in 
the two systems. Writing vector equation (C.l') in component 
form yields : 



v = v 

mx mx 

v 1 = v 

my my 

v 1 = v 

mz mz 

keeping in mind the definitions 
Figure Cl. Equations (C.2) can 
the direction cosines: 



- v ; (C.2a) 

P 

(C.2b) 

(C.2c) 

of the two systems given in 
be rewritten in terms of 



v 1 cos a' = v cos a - v ; 
m m p 



v' cos 6' = v cos 6 ; 

m m 



v' cos y 1 = v cosy, 
m ' m 1 



(C . 3a) 
(C.3b) 
(C.3c) 
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Squaring equations (C.3) and adding yields: 



22 2 2 2222 
(cos a' +cos S' + cos y') = v m (cos a +cos S +cos y) 



+ v - 2v v cos a 
P m P 



(C . 4 ) 



Since the sum of the squares of direction cosines is iden- 
tically equal to one, equation (C.4) becomes 



v 



.2 _ 



m 



2 2 

v + v - 2 v v cos a 
m p m p 



(c.5: 



From equation (C.3c) 



. ,cos y . 

v m = v' } ' 



m cos y 



(c.6: 



which, when substituted into equation (C.5), yields 



2 .cos y \ 2 

V ( —r) 

m cos y 



2 2 

v + v 
m p 



2 v v cos a 
m p 



(C . 7 ) 



Solving for | cos y * I ^ gives 



cos y 




cos 



v o 

a + (—2.) 2 



(C . 8 ) 



m 



Equation (C.8) is the fundamental result of this derivation. 
When substituted into equations (II. 5) and (II. 6), the result- 
ing expressions for P , equations (II. 7) and (II. 8), have 
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the advantage that they are written in terms of the frame 
of reference most often utilized when studying the rings. 
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APPENDIX D 



RELATING ORBITAL AZIMUTH TO THE DIRECTION COSINES 

A 

Consider a unit vector, u, as shown in Figure Dl . Its 
direction can be specified by the direction cosines (a,8/y)/ 
or by (y,9) where 0 is called the "orbital azimuth." How 
are these two representations related? 



z u 




Figure Dl . Angles Used to Specify the Direction of u 



First, note that 



2 2 2 
cos ex + cos 8 + cos y = 1 , 



or 



2 

cos 



e 



i 2 2 

1 - cos a - cos y . 



(D.l) 
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The x and y components of u are 



By definition 



u = cos a , 
x 



(D. 2a) 



u = cos 8 . 

y 



(D.2b) 



u 

tan 0 = — ^ , 

u 

x 



(D . 3 ) 



and substituting from equations (D.2) gives 



tan 0 = 



cos 8 
cos a 



(D . 4 ) 



Squaring equation (D.4) and substituting from equation (D.l) 
yields 



tan 2 6 = 1 - cos 2 y cos 2 Y 

cos a 



(D. 5) 



Rearranging equation (D.5) gives 



(1 



2 2 
+ tan 0)cos a 



1 



2 

cos y / 



but 



1 + tan 0 



2 

sec 0 



and 



1 - cos y 



sin y , 
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so 



sec 0 cos a = sin y 



Finally , 



cos a = sin y cos 0 , 



which is useful for expressing P c in terms of 



(D . 6 ) 



( D . 7 ) 

and 0 . 
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APPENDIX E 



SOME RESULTS OF CALCULATION OF P 

c 

This appendix consists of several examples of results 
for the calculation of P , based on the equation 




Calculations were made over the range of optical thickness 
pertinent to Saturn's A and B rings: 

0.1 £ x £ 2.0. 

The figures have been arranged to allow one to compare the 
effect of varying the parameters of the equation. 
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Figure El. Collision Probability for y = 0 
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Figure E2 . Collision Probability for y 
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Figure E3 . Collision Probability for y = 
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APPENDIX F 






ANNOTATED BIBLIOGRAPHY 

The purpose of this appendix is to briefly review those 
papers which should be of particular interest to one who 
wishes to study ring atmospheres. 



General Theory 

1. Blamont, J. 1974. "The 'Atmosphere' of the Rings of 
Saturn," in NASA SP-343 , edited by Palluconi and 
Pettengill, 125-129. 

Blamont 's paper is the transcript of his comments, 
delivered at a planning conference for the Pioneer II 
mission. His theoretical model is not detailed but 
results for the H 2 O, H, and OH atmospheres are given. 
Estimates of Lyman-a emissions are also presented. 

2. Dennefeld, M. 1974. "Theoretical Studies of an Atmos- 
phere around Saturn's Rings," in Exploration of the 
Planetary System, edited by Woszczvk and Iwaniszewska , 
471-481. 

Dennefeld was Blamont 's student. His paper predicts 
the density of the ring atmosphere (and that associated 
with Titan). The sources considered are sublimation, 
bombardment by meteoroids, solar and interstellar wind 
effects, and photodisssociation of H 2 O on the ring 
particles. The resulting estimates are outdated for the 
most part, primarily because the production mechanisms 
are now better understood . The conceptual model for 
calculation of atmospheric densities, however, is of 
enduring significance. It concisely gives a qualitative 
way to describe the atmosphere. One questionable asser- 
tion made in the paper is that the H atmosphere will 
be spherical due to collisions between the H atoms and 
ring particles; this is at odds with the results of 
the present thesis. Also, Dennefeld estimates the 
collision probability to be equal to the optical . 
thickness . 
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Ballistic Transport Models 



1. Durisen, R.H. 1984. "Transport Effects due to Particle 
Erosion Mechanisms," in Planetary Rings , edited by 
Greenberg and Brahic, 416-446. University of Arizona 
Press . 

Durisen presents a detailed model aimed primarily at 
prediction of "chip halo," i.e., macroscopic ejecta 
resulting from meteroid impact, formation (the model is 
in general applicable to the ring atmosphere, too). 

The assumption of "prompt absorption" (i.e., particle- 
molecule collisions are inelastic) is made for calcu- 
lations in his paper. In calculating collision 
probability, Durisen accounts for "slant path" of a 
molecule through the rings, but it is not clear whether 
he considered the relative speed of the particles versus 
the molecule. The results of his studies indicate that 
transport of eroded material across the rings may be 
significant for understanding the fine structure 
observed in the rings. 

2. Ip, W-H. 1983. "Collisional Interaction of Ring 

Particles: the Ballistic Transport Process," Icarus 

54: 253-262. 

Ip develops a transport model independently of Durisen. 
The results of the two different approaches are quite 
compatible. Ip is able to assume partially elastic 
collisions, since he employs the Monte Carlo method in 
his simulation. He takes the collision probability to 
depend on optical thickness only; a more detailed 
estimate of P c could easily be incorporated into his 
model, though. 



Sources (Production Mechanisms) 

1. Carlson, R.W. 1980. "Photosputtering of Ice and Hydrogen 
around Saturn's Rings," Mature 283: 461. 

Carlson argues that photodissociation of H 2 O on ring 
particles could account for the relatively high popu- 
lation of H atoms actually measured in the vicinity 
of the A and B rings. Note that he calculates collision 
probability using l-e _T with x equal to one. He takes 
the sticking coefficient to be 0.22, but this value is 
based on laboratory experiments (Brackmann and Fite, 1961) 
which do not closely match the conditions in Saturn's 
rings. The overall influence of the uncertainty in 
these two factors should not affect the order of magni- 
tude of his calculations, but probably limits the 
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precision of his estimation of total production rate 
needed to produce the observed H atom population. 

2. Cheng, A.F. and Lanzerotti, L.J. 1978. "Ice Sputtering 
by Radiation Belt Protons and the Rings of Saturn and 
Uranus," Journal of Geophysical Research 83: 2597-2602. 

This paper is a detailed summary of sputtering effects 
on ring particles . 

3. Cheng, A.F., Lanzerotti, L.J., and Pironelle, V. 1982. 
"Charged Particle Sputtering of Ice Surfaces in Saturn's 
Magnetosphere," Journal of Geophysical Research 87: 
4567-4570. 

Using information derived from the Pioneer II, Voyager I 
and Voyager II missions, the authors reconsider the 
sputtering mechanism, revising the result of their 
earlier (1978) paper. Since the actual ion fluxes 
observed were lower than those they assumed in 1978, 
the authors conclude that ion sputtering probably does 
not contribute significantly to the H atom atmosphere 
in the vicinity of the A and B rings of Saturn. 

4. Ip, W-H. 1978. "On the Lyman-alpha Emission from the 
Vicinity of Saturn's Rings," Astronomy and Astrophysics 
70: 435-437. 

This short paper suggests that H + ions escaping from 
Saturn's atmosphere may collide with ring particles 
where they are neutralized and re-emitted, thus 
contributing to the neutral H atmosphere in the vicinity 
of the rings. 

5. Morfill, G.E. et al. 1983. "Some Consequences of 
Meteroid Impacts on Saturn's Rings," Icarus 55: 439-447. 

Morfill et al. make careful estimates of the meteroid 
flux in the vicinity of Saturn. Using this estimate, 
coupled with a model for analyzing the effects of meteroid 
impact on ring particles, they proceed to make predictions 
about the ring halo, plamsa production, and the ring 
atmosphere. The authors conclude that meteroid impacts 
should contribute significantly to the ring atmosphere. 



Ring Characteristics 

1. Cuzzi, J. et al. 1984. "Saturn's Rings: Properties and 

Processes," in Planetary Rings , edited by Greenberg and 
Brahic, 75-199. University of Arizona Press. 

This paper is an up-to-date summary of information on 
Saturn's rings. It contains a short section devoted 
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to the ring atmosphere, which concisely explains 
current knowledge on that subject. Also, a special 
interest is an appendix which outlines radiative trans- 
fer theory (which is so important in deducing ring 
characteristics from observations) . 

2. Esposito, L.W. et al . 1984. "Saturn's Rings: Struc- 

ture, Dynamics and Particle Properties," in Saturn , 
edited by T. Gehrels and M. Matthews, 463-545. Univer- 
sity of Arizona Press. 

This paper nicely complements the preceding paper 
(Cuzzi et al. 1984). 

3. Marouf, E.A. et al . 1983. "Particle Size Distributions 

in Saturn's Rings from Voyager I Radio Occultation , " 

Icarus 54: 189-211. 

This detailed paper provides the first precise estimates 
of particle size distributions ever given for the rings 
(results are calculated for four locations in the rings) . 
The theory section in the paper also provides a clear 
application of radiative trasnfer theory. 

4. Weidenschilling , S. et al . 1984. "Ring Particles: 

Collisional Interaction and Physical Nature," in Planetary 
Rings , edited by Greenberg and Brahic , 367-415. Univer- 
sity of Arizona Press. 

Besides providing an in-depth analysis of the effect 
of collisions between ring particles on the evolution 
of the rings, this paper gives information useful in 
the study of particle properties. 



Miscellaneous 



1. Chandrasekhar, S. 1960. Radiative Transfer . Dover 
Publications, Inc. 

Chandrasekhar produced the definitive work in the 
field of radiative transfer (radiative transfer theory 
seeks to specify the radiation field produced when 
incident radiation is scattered and absorbed by media 
consisting of discrete particles). This book is rigorous 
in approach, in fact it can be appreciated as a work in 
mathematical physics, alone. It is not light reading. 

2. Pollack, J.B. 1975. "The Rings of Saturn," Space 
Science Review 18: 3-93. 

Pollack's paper is a review of knowledge on Saturn's 
rings up to 1975. Much has been discovered since then, 
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rendering invalid some of the predictions in the paper. 
This does not distract from the value of his work as 
a historical review. Reading this paper is an inspira- 
tional experience, because Pollack clearly demonstrates 
how knowledge about the rings increased through appli- 
cation of the scientific method. 

3. van de Hulst , H.C. 1981. Light Scattering by Small 
Particles . Dover Publications, Inc. 

This book is very readable. It essentially connects 
the mathematical physics of Chandrasekhar (1960) with 
the world of applied research. It provides the means 
by which the measurement of light scattering can be 
used to deduce the properties of the scattering medium. 
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